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Upon increasing significantly the nuclear elongation, the beta-decay energy grows. This paper 
investigates within a simple yet partly microscopic approach, the transition rate of the /!“ decay of 
the ^®^Cf nucleus on the way to scission from the exit point for a spontaneous fission process. A rather 
crude classical approximation is made for the corresponding damped collective motion assumed to 
be one dimensional. Given these assumptions, we only aim in this paper at providing the order 
of magnitudes of such a phenomenon. At each deformation the energy available for /3“ decay, is 
determined from such a dynamical treatment. Then, for a given elongation, transition rates for the 
allowed (Fermi) 0"^ —y O"*" beta decay are calculated from pair correlated wave functions obtained 
within a macroscopic-microscopic approach and then integrated over the time corresponding to the 
whole descent from exit to scission. The results are presented as a function of the damping factor 
(inverse of the characteristic damping time) in use in our classical dynamical approach. For instance, 
in the case of a descent time from the exit to the scission points of about 10“^° second, one finds 
a total rate of beta decay corresponding roughly to 20 events per year and per milligram of ^®^Cf. 
The inclusion of pairing correlations does not affect much these results. 


tio p = R 12 /R 0 of the distance i?i 2 between the mass 
centre of nascent fragments and the radius Rq of the cor¬ 
responding sphericai nucieus). One sees cieariy that the 
beta stabiiity vaiiey shifts from Caiifornium aimost to 
Nobeiium when the nuciear shape varies from a sphere 
to configurations ciose to the scission point where one 
has typicaiiy in such very heavy nuciei p ~ 2.5. 

This resuit couid be aiso attained upon considering the 
opposite behavior of proton and neutron separation en¬ 
ergies as couid be seen from an extremeiy crude iiquid 
drop approach. Indeed, let us consider a Liquid Drop 
energy formula comprising merely volume, surface, di¬ 
rect Coulomb and volume symmetry terms. For a spher¬ 
ical shape, one gets for the (positive) separation energies 
for protons {Sp) and neutrons [Sn) with a usual nota¬ 
tion and approximating the energy differences as differ¬ 
entiable quantities (i.e. making a continuous approxima¬ 
tion for the variables Z or TV): 

Sp = -\- ioc Z‘^A~^ — |as A~i — 2ac ZA~i 

Sn = -j- \ac Z'^A-^ - |as A-i 

( 1 ) 

where Sp^'^ and Sn'^ are deformation independent terms 
given e.g. for the protons by 

s(o, = p) 

It appears clearly that upon deforming the nuclear 
shape, both the surface and Coulomb deformation- 
dependent terms present in Sn will contribute to its de¬ 
crease. These terms are also present in Sp but are supple¬ 
mented by a term which increases with the deformation. 
We thus conclude that the difference Sp — Sn is increas¬ 
ing with deformation from about zero for the ground 
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I. PURPOSE 

Due to the decrease of the Coulomb energy and the 
increase of the surface energy, upon strongly deform¬ 
ing a nucleus, it may happen that a nucleus which 
is /3“-stable at equilibrium deformation gets a positive 
Q{/3~) = {M{Z,A) — M{Z -b 1, A) — me)c^ at large de¬ 
formation. This is exemplified on Fig. [T] within a Liq¬ 
uid Drop Model (LDM) approach [I|, where the masses 
of A = 252 isobars (from Plutonium to Kurchatovium) 
are plotted for various deformations (defined as the ra- 



FIG. 1. Mass excess of Pu - Ku nuclei with A=252 evaluated 
within the LSD model [J]. Each curve corresponds to different 
elongation of nuclei p defined in the text. These values are 
reported on the r.h.s. side of the every second curve. The 
upper curve (labelled as ”sph.”) presents the mass excesses 
calculated for spherical nuclei. 
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FIG. 2. Proton (p) and neutron (n) single-particle levels in 
the ground-state (l.h.s.) and at the scission point (r.h.s.). 



FIG. 3. Macroscopic (solid line) and macroscopic- 

microscopic (dashed line) fission barriers and the beta-decay 
energy Qjs (dotted line) of ^®^Gf as a function of nuclear elon¬ 
gation. 


state of a /3 stable nucleus to a value which could be 
of the order of some MeV. Resorting now to state of the 
art macroscopic-microscopic calculations (using the LDM 
model of Ref. [ij and the Yukawa folded potential of 
Ref. i) including pairing correlations using a seniority 
pairing force, one gets an inversion of the relative posi¬ 
tion of neutron and proton Fermi levels of the ^®^Cf beta 
stable nucleus upon increasing the elongation, the lat¬ 
ter becoming larger than the former. This is illustrated 
on Fig. [2] showing the neutron and proton single particle 
(s.p.) spectra around the Fermi levels for p r; 0.9 and 
p R! 2.8 corresponding respectively to the ground state 
and the scission point solutions. 

We are therefore considering the following /3“ decay 
during the spontaneous fission process (beyond the exit 
point): 

2"2cfi54 Esi53 + e-+D (3) 

We limit ourselves to allowed transitions. The spin and 
parity of the compound nucleus being conserved during 
the fission process (in the absence of particle emission) we 
are considering thus a O'*' —O'*" transition. This rules 
out a Gamow-Teller transition and we are thus simply 
left with a Fermi transition. 


II. APPROXIMATIONS 

At this very preliminary stage of our study, we make 
the following approximations pertaining to the static and 
dynamical calculations as well as to the evaluation of 
transition rates. 

a) We consider only a single fission path assumed to be 
the most probable one out of static macroscopic - micro¬ 
scopic calculations of the above discussed type. Pairing 
is included, as it should, to define adequately this opti¬ 
mal path. On Fig. [31 we have plotted the corresponding 


fission barrier down to the scission point, along with the 
underlying LDM one. 

As a result we find that apart from the ground state 
region, and in particular at large deformations (namely 
R 12 /R 0 > 1-5), the LDM energy curve represents a rea¬ 
sonable approximation of the results, where shell and 
pairing effects are taken into account. However, for 
the calculated ground state one has definitely to make 
a choice about its exact definition. One may rule out 
the LDM solution since shell and pairing effects increase 
its height. As a good approximation of the macroscopic- 
microscopic exit point, as can be checked on Fig. [31 is the 
exit point on the LDM fission barrier (on the side of the 
descent towards fission) which is located at exactly 4.8 
MeV below the saddle point. This value corresponds to 
the experimental fission barrier height (noting en passant 
that our macroscopic-microscopic calculations do repro¬ 
duce it quite well). 

In what follows we will use also, instead of the ratio 
p = R 12 /R 0 , a collective variable q obtained by a trans¬ 
lation of p such that g = 0 at the exit point. We define in 
a somewhat arbitrary fashion a single scission point. It 
corresponds to the point on the LDM fission path, where 
the nuclear liquid drop ceases to be consisting of one sin¬ 
gle piece (it is obtained here for p r; 2.9). It corresponds 
to a gain in energy of AEmax ~ 28 MeV with respect to 
the exit point energy (or the ground state one, since they 
are of course the same). 

On Fig.|3]we have also plotted the LDM Q{f3~) value 
as a function of p. It varies from -1.2 MeV at sphericity 
to enter the instability region for p r: 1.5 which happens 
coincidentally, to be very close to the exit point to reach 
r; 2.5 MeV around the scission point. 

b) In the present approach, which aims only at provid¬ 
ing orders of magnitude, we make the following approxi¬ 
mations for the collective dynamics: 

- As a function of the collective coordinate q, one may 
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approximate reasonably (see Fig. IS the potential energy 
Ep{q) as a linearly decreasing function when q varies from 
zero to qs which the value of this parameter at scission. 
One has thus 

Epiq) = -aq a Ri 20 (0 < g < (4) 

- We approximate the mass parameter m(g) as being con¬ 
stant m(q) = fjL during the descent from the exit point to 
scission and equal to the reduced mass between the two 
final fragments (to assess the relevance of this approxi¬ 
mation see Ref. i). 

- We assume that the damping of the motion results from 
a simple friction force —kq with k > 0. 

Then the equation of motion yields the following triv¬ 
ial solution (introducing of course a vanishing collective 
velocity at the exit point) 

Qf —Mm /«\ 

")], (5) 

which could be numerically inverted to get the time t{q) 
as a function of the deformation q. 

Then the collective kinetic energy will be given by 

EK{q) = ^ . (6) 

The available energy is 

^E{q) = -Ep{q) = aq , (7) 

from which one gets the excitation energy E*{q) as 

E*{q) = ^[x-g{x)-^g{xr] , ( 8 ) 

where the dimensionless quantity x is defined by 


and with 


III. EFFECTIVE Qp VALUES AND AVAILABLE 
PHASE-SPACE 


At a given elongation q or at given time t(g), the parent 
nucleus has an excitation energy of E*{q). The effective 
Qp{q) value for a transition from a neutron single particle 
state i to a proton single particle state / is thus given by 


Qf^\q) = E*{q) + en{i) - epif)) + 0.78 MeV , (12) 

where the last constant origins from the mass difference of 
neutron and proton plus the electron mass. Defining the 
neutron Fermi level Xn and similarly for protons W, one 
will consider an effective neutron Fermi level as A„ (g) as 


Af (g) = A„ + E*{q) + 0.78 MeV . 


(13) 


Possible transitions between neutron s.p. states and pro¬ 
ton s.p. states will require the following conditions to be 
satished 


(*) = Kpif) , 

>p < ep(/) < e„(*) < Xf{q) 


IV. TRANSITIONS RATES 

The nuclear matrix element for such a {i, f) Fermi 
transition is given with a usual notation by 

Mij = JJJ . (15) 

A special note is to be made here, about the intrinsic 
parity breaking. Let us call Idt) the intrinsic wavefunction 
breaking the left-right symmetry, either for the parent or 
the daughter states. One should project a positive parity 
state out of it, as 

|vi/(+))| = 4(l'i')l + n|4^)), ( 16 ) 


g{x) = 1 - e “ . (10) 

c) To compute beta decay rates, the effect of pairing 
correlations will be taken into account. For the sake of 
clarity we will first discuss our approach without them 
and in a separate section the modifications brought in by 
their inclusion. 

The conservation of the axial symmetry is assumed 
along the path. Each s.p. state is thus defined in par¬ 
ticular by the usual angular momentum projection K 
quantum number. 

We will assume that there is no polarization i.e. that 
the s.p. wavefunctions in the ^®^Cf or ^®^Es mean helds 
are the same and furthermore that one has, with a trans¬ 
parent notation, the following relation between the bind¬ 
ing energies of the parent and daughter nuclei (Koop- 
mans approximation) 

E(252Es) = F;(252cf) -H Cp - e„ . (11) 


where 11 is the parity operator. 

The above nuclear matrix element will thus comprise 
two parts corresponding to the two overlaps ('I'(i)|d>(/)) 
and (4t(i)|n|'I'(/)). However it is known from parity pro¬ 
jection calculations of the fission barrier of heavy nuclei, 
see Ref. Q, that somewhere before the second fission 
barrier and beyond, the intrinsic parity breaking defor¬ 
mation is so large that the overlap between the corre¬ 
sponding wavefunction and its parity image is negligible. 
We are just left here with the ('I'(i)|tI/(/)) overlap as in 
the non-parity breaking case. Yet, in this case one can¬ 
not, of course, apply a selection rule on the quantum 
number tt which is not conserved but only on the angu¬ 
lar momentum projection K which we have assumed to 
be a good quantum number. 

The associated transition rate is @ 


i?,./ = G" 


mlc^ 

2tt^ h7 


fiZ = 99 ,Qy) , 


(17) 
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which after inserting of the week interaction coupling 
constant (G) and other constants values takes the fol¬ 
lowing form; 

= 1.105 10-4 |M,j|2 f{Z = 99,Qy). (18) 

A rough estimate of the function / as a function of the 
effective Q value Qp’^\q) is given for Z = 99 by 

logio(/) = 3-5 logio(Q^*’^^(g)) +3 > (19) 

where (q) is expressed in MeV. 


V. INCLUDING PAIR CORRELATED 
NUCLEAR STATES 


When pairing correlations are included in both the par¬ 
ent and daughter nuclear states (assuming again no po¬ 
larization effects, i.e. taking the s.p. wavefunctions, oc¬ 
cupation probabilities and quasi-particle energies as ob¬ 
tained in the parent nucleus) one has just to multiply 
the overlaps Mij by the BCS factor Upi^f)Vn(i) while the 
effective {q) becomes 


Qf^\q) = E*{q) - [Ar(*) + A^n/)] + (A„ - Xp) 
-fO.78 MeV , 


( 20 ) 

where E‘^(i) and E^P{f) are the initial neutron and the 
final proton quasiparticle energies. 

In this case, one may consider all transitions from sin¬ 
gle neutron to single proton states which satisfy the fol¬ 
lowing conditions beyond the angular momentum com¬ 
ponent matching rule (Kn(i) = Kp{f)): 

- the proton state should not be fully occupied (i.e. lying 
below the valence space), 

- the neutron state should not be fully unoccupied (i.e. 
lying above the valence space), 

- the final state energy should be lower or equal to the 
initial state energy, i.e. {q) > 0 . 


VI. TOTAL TRANSITION RATES 

At a given time (or at a given elongation) the total 
rate of decay will be obtained by summing all individual 
rates Rij to get R{q). For the whole descent from exit 
to scission, one will get the number of decay for a sin¬ 
gle fission, i.e. the probability of decay, by integrating 
R{q) over the whole time of descent from the exit to the 
scission points 

Pdecay = f R{q) -{l-e~E) dq. (21) 
Jo a 

This probability, and similarly the pre-scission kinetic 
energy E^ = Ex^qs), are of course strongly dependent 
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FIG. 4. Collective elongation q = {R 12 — R’u^) / Ro as function 
of time (eq.[5]) for different values of the nuclear damping k/fi. 


on the retained value for the friction parameter k. In the 
following we will use our estimates of Pdecay as a function 
of the time tsc to evaluate the number of /3-decays from 
a sample of ^^^Cf in a given time period. 


VII. RESULTS 


t [10“^^ s] 
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FIG. 5. Pre-scission kinetic energy for a few values of the re¬ 
duced friction as function of time. The time at which scission 
conhguration is reached for a given fc/p value is marked by a 
cross. 

Let us first use the results of our simple dynamical 
calculations to describe the collective motion from the 
exit to the scission points. Our results will be presented 
for a few values of the reduced friction parameter 
which inverse is the damping time ©• The damping 
parameters fc//i is given in units MeV//i ~ 1.5 

As seen in Fig. 01 the collective descent from g = 0 to 
q Ki 1.4 takes from Ri 3.5 IQ-^^g for a large dumping time 
(fc//r = 1 in the above discussed units) to ~ 18 IQ-^^g for 








5 


a small damping time {k/fi = 6)). In all considered cases 
the collective motion is completely damped at scission 
as exemplified on Fig. [5l The pre-scission kinetic energy 
Ek values range from a couple hundred keV for the large 
damping case (fc//r = 6)to~5 MeV for the low damping 
case (k/^ = 1). 

The resulting number of beta decays for a 1 mg sample 
of is given in Fig. [51 as a function of the descent 

time tsc from the exit to the scission points (or equiva¬ 
lently as a function of the damping parameter k/fi). It 
ranges from 6 for the low damping case {k/ji = 1) to 
« 46 for the large damping case {k/fj, = 6). The range of 
considered values of the damping parameter fc//i is taken 
following estimates made in Ref. Q (see Fig. 10 of this 
Reference). 


k/p [MeV/h] 

0 1 2 3 4 5 6 



FIG. 6. Number of beta decays per year for a 1 mg sample 
of as a function of the descent time tsc from the exit to 

the scission points. 

The previous figures are given for an evaluation which 
does not take into account pairing correlations effects. 
As seen of Fig. IH including these correlations increases 
slightly the above given numbers (only by ~ 2% for 
fc//i = 6), this correction being an increasing function 
of the damping parameter fc//i. 


VIII. CONCLUSIONS 

This paper aimed at pointing out the strong depen¬ 
dence of the beta-decay stability as a function of the 
nuclear deformation. From state of the art macroscopic- 
microscopic calculations it has been claimed that a very 
heavy nucleus like while beta-stable as well-known 

in its ground state becomes instable near the fission exit 
point. Through a very simple description of the collective 
dynamics, the transition rates of the Fermi beta decay 
during the spontaneous fission process of this nucleus up 
to the scission point, have been calculated. They depend 
of course on the amount of damping of the fission collec¬ 
tive mode which regulates the excitation energy available 
at each time, for such a decay. In so far as these rates 
could prove to be experimentally reachable, they would 
provide a much needed source of informations on the pre¬ 
scission dynamics, as e.g. the average pre-scission kinetic 
energy or the value of the reduced friction parameter. 
Taking or not into account pairing correlations is not an 
important factor. 

As a result a relatively weak number of beta decay 
are expected from our results during the descent from 
the exit to the scission points for standard values of the 
damping parameter. The capacity of existing experimen¬ 
tal devices to assess this conclusion for tractable amounts 
of ^®^Cf is questionable. Improving these so far very 
crude estimates by a better treatment of the static and 
dynamic parts of these calculations as well as exploring 
this phenomenon for other nuclei (involving possibly the 
inclusion of Gamow-Teller decays as well) will be under¬ 
taken. 
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